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Abstract. The bottleneck network flow problem (BNFP) is a generalization of several well-studied 
bottleneck problems such as the bottleneck transportation problem (BTP), bottleneck assignment 
^>-^ . problem (BAP), bottleneck path problem (BPP), and so on. In this paper we provide a review 

' of important results on this topic and its various special cases. We observe that the BNFP can 

, be solved as a sequence of O(logn) maximum flow problems. However, special augmenting path 

04 ' based algorithms for the maximum flow problem can be modified to obtain algorithms for the 

^ ' BNFP with the property that these variations and the corresponding maximum flow algorithms 

^ , have identical worst case time complexity. On unit capacity network we show that BNFP can be 

' solved in 0(min{m(n log n)^ , m' ^log n}). This improves the best available algorithm by a factor 

of \/\og n. On unit capacity simple graphs, we show that BNFP can be solved in 0(m\Jn log n) 



time. As a consequence we have an 0(rny/n log n) algorithm for the BTP with unit arc capacities. 
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. 1. Introduction 

. Let E) be a directed graph such that |y| = n and \E\ = m. For each arc G E, a weight 

00 I Cij and a capacity Uij are prescribed. Also, for each node i £ V an integer 6j is associated. Then 

■ the bottleneck network Bow problem (BNFP) can be formulated mathematically as follows: 



BNFP: Minimize maxjcjj : Xij > 0} 

Subject to 

i - Yl. = bi yi £V 

O ■ {j-{i,j)&E} {j:{j,i)eE} 

< Xij < Uij V(i,i) G E 

^ I Here Xij is the flow on arc We assume that Cij and Uij are integers for all G E and there 

■ are no multiple arcs in G. The integer number bi associated with node i represents the supply or 

demand at the node. A node i is called a supply node if ftj > and a demand node if hi < 0. If 

hi = 0, then i is a transshipment node. We assume that ^ ftj = 0. 

i&V 

An interpretation of the bottleneck network flow problem can be given as follows: Suppose that 
perishable goods are shipped from the supply nodes to the demand nodes. The goods can be stored 
without damage at the nodes but are perishable in time on transit. Assume that Xij is the amount 
of goods shipped along the arc and Cij is the shipment time along the arc Then the 

BNFP objective function measures the largest time a shipment in transit and we want to minimize 
this time. 
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To the best of our knowledge, the general form of the bottleneck network flow problem has not 
been studied in literature except the a generalization on algebraic flows [6]. However, many special 
cases of this problem are well-studied. One of the most well studied special cases of BNFP is the 
bottleneck transportation problem (BTP). In this case, the graph G is bipartite with the generic 
bipartition of V as F = Vi U V2, such that i £ Vi implies bi >0 and i G V2 implies bi < 0. Most of 
the literature on bottleneck transportation problem assumes that the arcs are without capacities 
(i.e. Wij =oo)[8l[IIl[ia[2ll[H[3a[Ml[Sl[Ml[371[Ml[121[^ although some papers 

admit finite capacities. 

In BTP, if 6i = 1 for i G Vi, bi = -1 for i £ V2, \Vi\ = IV2I and Uij = 1 for all (i, j) G E, then the 
BTP reduces to the well known bottleneck assignment problem (BAP) [HHl[5l [T0l[T2t [T7 1[T8| [T9l 

Another well studied special case of BNFP is the bottleneck path problem(BPP). Let s and 
t be two specified vertices in G. If 6(s) = l,b{t) = —l,bi = for i G 1/ \ {s,t} and Uij = 1 for 
every G E, the resulting bottleneck network flow problem is called bottleneck path problem 
(BPP) [IlEniliaiillllSllIlllIalllS]. Further, BNFP is a special case of the bottleneck linear 
programming problem (BLP) [22 ^ 1551 [62] . 

The paper is organized as follows. In section 2 we provide a literature survey on existing algo- 
rithms for various special cases of BNFP. In section 3 we discuss some basic algorithms for solving 
the BNFP. We first observe that BNFP can be solved as an O(logn) sequence of maximum fiow 
problems. We then identify special maximum flow algorithms that can be modifled to solve BNFP 
with the same worst complexity as that of solving just one maximum flow problem by these al- 
gorithms. Section 4 deals with the unit capacity networks. We flrst observe that the best known 
maximum flow algorithm for unit capacity networks can be extended to handle arbitrary capacities 
on arcs incident on source and sink nodes. We then show that BNFP on unit capacity graphs 
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can be solved in 0(min{(n • logn)3 • m, ma • yTogri}). This improves the best known algorithm 
for BNFP on unit capacity networks by a factor of ^log n. For unit capacity simple networks, we 
obtain a complexity bound 0{myjn log n). This algorithm can be viewed as a generalization of the 
algorithm of Gabow and Tarjan for the BAP [6OJ and also provides an 0{m^/n logn) time bound 
for the bottleneck transportation problem with unit capacities. 

2. Literature Review 

In this section we provide a review of known results on BTP, BAP and BPP which are special 
cases of BNFP. To the best of our knowledge no review papers on the topic is available. We keep the 
review section brief, highlighting only important results. For details, the reader is referred to the 
original papers. Let us first consider the BTP. Most of the known algorithms for this problem can 
generally be categorized into three groups: (1) primal algorithms (2) augmenting path algorithms 
and (3) threshold algorithm. Primal algorithms start with a feasible solution and try to find a better 
solution. Since the different objective function values of BTP solutions are at most m = \ Vi\ x IV2I, 
the number of improvement steps is 0{m). Algorithms discussed by Hammer |331l34j. Garfinkel and 
Rao [21], Bhatia, Swaroop and Puri [38] etc. falls in this category. Another class of algorithms build 
a solution by means of augmenting paths. Algorithm proposed by Derigs and Zimmermann [T3] 
is an example of such an algorithm. The algorithm augments fiows along bottleneck paths until a 
feasible(and hence optimal) solution is obtained. The complexity of this algorithm can be verified 
to be 0{S ■ f{m,n)) where S = J2ieVi 1^*1 fi^,n) is the complexity of BPP. The threshold 
algorithm sets a threshold for the optimal objective function value and checks the existence of a 
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solution satisfying this threshold value. Depending on the outcome, the threshold value is adjusted 
and the process is continued. Some of the primal algorithms can also be viewed as a threshold 
algorithm. For a discussion on threshold algorithms for general combinatorial bottleneck problems 
we refer to the paper by Edmonds and Fulkerson [15]. When the number of supply (demand) nodes 
are fixed, say k, Hochbaum and Woeginger |36] showed that the BTP can be solved in 0{n) time. 
A Special case of this problem when k = 2 has been studied by Varadarajan [61] who gave an 0{n) 
algorithm and Szwarc [59jwho proposed an 0{n ■ logn) algorithm. Many of the works on BTP are 
relatively old and these papers do not discuss complexity results. It is easy to obtain a binary 
search version of the threshold algorithm to solve BTP as a sequence of O(logn) maximum flows 
in a bipartite graph. This result extends in a straightforward way to obtain a threshold algorithm 
for the BNFP which also solves O(logn) maximum flows. 

Most of the algorithms known for the bottleneck assignment problem (BAP) can also be cate- 
gorized as primal algorithms [31j, augmenting path algorithms [12] and threshold algorithms [23] . 
These algorithms can be viewed as specializations of the corresponding algorithms for BTP. The 
best known algorithm for BAP is a hybrid algorithm that uses a binary search based threshold 
algorithm together with an augmenting path algorithm. Using this approach, Garbow and Tar- 
jan [60] obtained an algorithm of complexity 0{m\Jn logn) to solve BAP. In the threshold phase 
of this algorithm, a "relaxed" feasibility problem is considered to obtain a partial solution, which 
is extended into a full solution by means of augmenting paths. It is the best known time bound 
for BAP on sparse graphs. Using a similar approach, Punnen and Nair [50] proposed an 0{n^/nm) 
algorithm by considering a slightly different "relaxed" problem. This bound is the best known 
for solving BAP on dense graphs. When the arc weights are Euclidian distances, Efrat, Itai and 
Katz [1] proposed an 0{m) algorithm. When V C M'^, Efrat and Katz [17] proposed an 0{n^'^) 
time algorithm for d < 6, and a subquadratic time algorithm for d > 6. If the underlying norm 
is Loo, then the complexity bound is O(n^'^log^'^n) for d > 2. When Cij = ai ■ bj, Eiselt and Ger- 
chak [18j proposed a simple non-iterative scheme. Probabilistic results on BAP are discussed by 
Pferschy [47] and specially structured cost matrices are considered by Cechlarova [10], Eiselt and 
Gerchak [18]. 

A natural approach to solve the bottleneck path problem (BPP) is to consider modifications of 
the shortest path algorithms. Many authors considered modifications of the Dijkstra's algorithm 
for shortest path [14j to solve BPP [51]. The complexity of such an algorithm is O(n^) for a 
straightforward implementation. Using Fibonacci heaps, the method can be implemented in 0{m + 
nlogn) time [3]. Listrovio and Khrin [H] also proposed a related algorithm explained using s — t 
cuts. Their algorithm starts from an s — t cut K = [S, S], where S = {s}, S = V — {s}. The 
maximum capacity of this cut is set to be a lower bound of the objective function value, and K is 
iteratively modified by increasing S and decreasing S until the sink node t € S. Fernandez, Garfinkel 
and Arbiol [52] presented a binary search based threshold algorithm [52] for BPP. This paper also 
discusses an application of BPP in the context of combining (mosaicking) two or more aerial 
photographs into a single image in the production of photographic maps. Inspired by an algorithm 
of Gabow and Tarjan for bottleneck arborescence problem, Punnen [49] showed that if a bottleneck 
combinatorial optimization problem of size m with ordered weights can be solved in 0(,^(m)) time, 
then the problem with arbitrary weights can be solved in 0(^(m) log*(m)) time, where log* n is the 
iterated logarithm of m. As a consequence, the BPP can be solved in 0{m log* m) time. Combining 
this approach with modification of Dijkstra's algorithm discussed earlier, which uses Fibonacci 
heaps, the best known complexity for BPP on a directed graph is 0(min{m + nlogn, m log* m}). 
Georgiadis [25] showed that BPP can be solved in 0{T{m)) time where T{m) is the time for sorting 
the edge costs of the underlying graph. BPP on an undirected graph is simpler and can be solved 
in linear time using a binary search based threshold algorithm coupled with subgraph contractions. 
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For details of this algorithm, we refer to Punnen [48j . Sensitivity analysis for BPP have been 
investigated by Ramaswamy, Orlin and Chakravarty [53] . It is easy to show that all pair bottleneck 
path problem on an undirected graph can be obtained by computing just one minimum spanning 
tree. 

3. Basic algorithms for BNFP 

The BNFP can be formulated as a minimum cost flow problem with exponentially large arc 
costs j32] . However, this is not a practical approach to solve the problem. We now consider some 
basic algorithms to solve BNFP, which are generalizations of the corresponding algorithms for the 
bottleneck transportation problem (BTP). 

For any real number a, let G{a) = {Va,Ea) denote the spanning subgraph of G with Va = V 
and Ea = £ E : Cij < a}. Let Co-(i) < Co-(2) < • • • < Co-(<^) be an ascending arrangement of all 

distinct arc weights of G. Let 6 = J2 — where S is the collection of supply nodes and T 

is the collection of demand nodes. 

The auxiliary graph G* = {V*,E*) corresponding to any graph G{V, E) is defined asV* = VU {s, t}, 
E* = E U {{s, i) : i £ S} U {{j, t) : j £ T}, and s ^ V , t ^ V . Here s is called a source node and t 
is called a sink node in G*. The capacity Usi of arc (s, i) is bi for all i £ S and the capacity Ujt of 
arc (j, t) is —bj for all j £ T. The weights Csi of arc (s, i) for all i E 5" and Cjt of {j, t) for all j £ T 
are set to be 0. Clearly BNFP is feasible if and only if G* has an s — t flow of value 6. 

Theorem 3.1. Assume that the BNFP is feasible. Let k £ {1, 2, . . . ,ip} be the smallest index such 
that the maximum flow in G*(co-(yt)) is 6. Then any flow x^ in G*{c„(^i.-j) of value 6 provides an 
optimal solution x to BNFP. 

The straightforward proof of Theorem 13.11 is omitted. Note that the solution x in the above 
theorem is obtained by simply dropping the flow values on arcs incident on s and t from x". 

Based on Theorem l3.1l we see that the if we get the value of k, then we can get an optimal solution 
for BNFP. In fact, k can be identifled by using different search strategies. Using binary search over 
the set {co-(i), Co-(2)) • • • 1 Co-(ip)}) by Theorem 13.11 it can be verifled that BNFP can be solved by 
solving 0(logip) = O(logn) maximum flows. We call this algorithm the binary search threshold 
algorithm. This observation raises an interesting question: "Is it possible to solve BNFP using less 
than O(logn) maximum flow computations?" We do not have an answer to this. Later in section 
4 we will see that for unit capacity networks, we can solve BNFP using 0(\/logn) approximate 
maximum flow computations. For the general BNFP, let us consider a closely related question: "Is 
it possible to modify a maximum flow algorithm to solve the BNFP within the same time bound as 
that of solving the maximum flow problem with the original algorithm?" As we show below, this 
is doable in some cases and the question is open for other cases. 

Perhaps, the simplest such example is the generic augmenting path algorithm for maximum 
flows [3j. Here, we start with the graph G*(co-(i)) and augment flows from s to t by augmenting 
paths. If a flow value of 5 is reached, we have an optimal solution to BNFP. Otherwise we add 
to G*(co-(i)) all arcs of weight Co-(2) to obtain the graph G*(co-(2)) and search for augmenting path 
is continued. Continuing this process by adding new classes of arcs in the increasing order of 
weights and the algorithm terminates when a flow x^ of value 5 is identifled. By Theorem 3.1, an 
optimal solution to BNFP can be recovered from x^. As in the case of the generic augmenting 
path algorithm for maximum flows, the complexity of this algorithm is 0{m6 + (plogip) = 0{mnB) 
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where B = max IbA. This algorithm is a variation of the augmenting path algorithm of Derigs 

and Zimmerman [13j designed for BTP. One major difference is that we do not use bottleneck 
path computations to identify augmenting paths which results in slightly improved complexity 
for this more general problem. Further the algorithm provides a natural linkage with the generic 
augmenting path algorithm for maximum flows. We call this algorithm the BNFP augmenting path 
algorithm. 

Let us now consider the maximum capacity augmenting path algorithm for maximum flows [TB] . 
This algorithm has polynomial complexity and we show that the algorithm can be easily modified 
to solve BNFP without increasing the complexity bound. Rather than considering maximum ca- 
pacity augmentations, we augment flows along paths with large enough residual capacity to avoid 
maximum capacity path computations. This is possible because our target the maximum flow value 

5 is known a priori. We call the resulting algorithm for BNFP the large capacity augmenting path 
algorithm. 

Let G = iy , E) be a subgraph of G* containing s and t and Gr{x) be the residual graph [3] with 
respect to an s — t flow x in G. 

Lemma 3.2. // G has a maximum flow from s to t of value 5, then for any flow x in G, there 
exists an augmenting path in Gr{x) with residual capacity at least where m* = \E*\ and 

v{x) is the value of flow x. 

Proof. Since Gr{x) is the residual graph of G with respect to flow x, then in Gr{x), a flow of value 

6 — v{x) can be represented as path flows. By the Flows Decomposition Theorem [3] {pp 79-81), at 
most m s-t paths have non-zero flow, therefore there must be an augmenting path in Gr{x) with 

residual capacity at least ^— Since m* > m, ^ > ^ and the result follows. □ 

Lemma 3.2 implies that we can modify the augmenting path algorithm by performing each 
augmentation along a path with capacity at least — 7^1-^- If no such path exits, we can safely 
conclude that G does not have a maximum s — t flow of value 5. Based on this idea, we present 
the large capacity augmenting path algorithm below: 

Algorithm large capacity augmenting path 
begin 

construct G* from G and solve a maximum flow problem on G*; 
if the maximum flow value < 6 then stop. BNFP is infeasible; 
else do 

let Co-(i) < Co-(2) < • • • < be an ascending arrangement of all 

distinct arc weights of G; 
k = 0, x^ = zero flow, v{x^) = 0; 
repeat 

k = k + l; 

let x'' := x^~^, v{x^) := v{x''~^); 
begin 

obtain the residual graph G*(co-(fc)) with respect to the 

flow x'^~^ from G*{c„(^i._i-^) by adding the arcs 

whose weight Cij = c^i^^-^ ; 
while (G*(co-(-fc)) contains an augmenting path P of 

residual capacity at least — ^) do 
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let e be the residual capacity of P; 
augment e units of flow along P; 
v{x'^) = v{x^) + e; 
update and G*(co-(fe)); 
end while; 
end if; 
until {v{x^) = 5) 

compute and output the optimal BNFP solution corresponding to x^; 
end; 
end; 

To test the conditions of the while loop of the above algorithm, we can construct a graph G'*(co-(A;)) 
by removing all the arcs in G*.{c„(^j.-^) whose residual capacity is less than ^ — -. Then G'*(co-(fc)) 
has an augmenting path of capacity at least ^ — - if and only if G'*(Co-(fc)) has an s — t path. 

To establish the complexity of the large capacity augmenting path algorithm, we prove the 
following theorem, which is a variation of a result by Edmonds and Karp [16] and Goldfarb and 
Chen [28j for the maximum capacity augmenting path algorithm for the maximum flow problem. 

Theorem 3.3. In the large capacity augmenting path algorithm, the number of augmentations is 
0(mlog5) = 0{mlog{nB)), where B = max|6j|. 

Proof. Let v{x^), f (x^), • • • , v{x^) = 5 be a sequence of flow values generated by the large capacity 
augmenting path algorithm. Thus k is the total number of augmentations performed. Assume 
di = v{x^~^^) — v{x'^) and = 5 — ^(x*), then 

di = v{x'+^) - v{x') = Ai - A,+i (1) 
Since the augmenting paths have capacity at least ^ , we have 



and hence Aj < m*di (2) 
From (1) and (2), 

A,+i< A,(l-;i^). 

Therefore 

Ap < Ai(l - < 5{1 - < ,5e-P/-* 
We want the largest p such that 

Ap > 1 . Thus p < m log 5 
Since v{x^) = 6, we have k < p + 1. Therefore k = 0{m* log 5) = 0{mlog{nB)). □ 

Theorem 3.4. The large capacity augmenting path algorithm correctly solves BNFP in 0{m? log{nB)) 
time. 

Proof. Starting with a zero flow in G*(co-(i)), the algorithm looks for the smallest index k such that 
G*(cCT(fc)) has a flow of value 6. At a typical iteration we have a flow x^ in G*(ca-(r)) for some r with 
value v{x'^) and v{x^) < 6. Then introduce arcs Co-(r+i) to obtain the graph G*(Co-(r+i)). Clearly x^ 
is a feasible flow in G*(co-(r+i))- If the residual graph of G*(co-(r+i)) does not have an augmenting 
path of value at least — ^r-^, by Lemma 3.2, G*(co-(j._|_i)) does not contain a flow of value 5 and the 
arc Co.(j,_|_2) is added to G*{c^(^r+i)) to obtain G* (00.(^+2)) and c^(^r+2) becomes a new lower bound 
for the optimal objective function value of BNFP. Otherwise the flow is improved by using an 
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augmenting path of capacity at least — . By Theorem 3.3 the number of augmentations is 
bounded by 0{mlog{nB)). The complexity of performing the augmentation step is 0{m), so the 
overall complexity of the algorithm is 0{m'^ log{nB)) . □ 

It may be noted that we could not obtain a variation of the shortest augmenting path algorithm 
for maximum flows that solves BNFP with in the same time bound as the corresponding maximum 
flow algorithm. 

4. BNFP IN Unit Capacity Networks 

On a unit capacity graph, it is well known that the maximum flow problem can be solved in 
0(min{m^/^, n^/^m}) time [20j . Suppose G is a unit capacity graph on which a BNFP is defined. 
Then the corresponding auxiliary graph G* = {V*,E*) will be of unit capacity except for the arcs 
incident on the source node s and sink node t. We first observe that the maximum flow problem 
in such a graph can also be solved in 0(min{m^/^, n^/^m}) time. 

A graph G with a source node s and a sink node t is called an almost unit capacity graph if (1) 
arcs incident on s or t or both have arbitrary capacities (2) all other arcs are of unit capacity and 
(3) any s — t path in G contains at least one arc which does not incident on s or t. The maximum 
flow problem restricted to an almost unit capacity graph is called almost unit capacity maximum 
flow problem (AMFP). An almost unit capacity graph G is simple if every node in G has at most 
one incoming arc or at most one outgoing arc. The corresponding maximum flow problem is called 
almost unit capacity simple maximum flow problem (ASMFP). 

4.1. Flows in Almost Unit Capacity Graphs. Let us now discuss the maximum flow problem 
in almost unit capacity graphs and almost unit capacity simple graphs. Our algorithms are similar 
to the unit capacity maximum flow algorithm of Edmonds and Karp [16] as discussed in Ahuja and 
Orlin [3J. We only discuss the primary results without proof. An interested reader could construct 
the proofs with appropriate modifications of the corresponding unit capacity case or can find it in 
the thesis [63] where details of the algorithms of this paper are given. 

It is well known that the shortest augmenting path algorithm [3j solves the maximum flow 
problem in 0{mn'^) time. On unit capacity graphs, the complexity can be reduced to 0{mn) [3]. 
The shortest augmenting path algorithm maintains distance labels that are non-decreasing and 
terminates when the distance label d{s) of node s satisfies d{s) > n. Let D < nhe a parameter. In 
the shortest augmenting path algorithm, if we discard all nodes i with distance label d{i) > D from 
further consideration, we get an approximate version of the shortest augmenting path algorithm. 
We refer to this algorithm the D-shortest augmenting path algorithm. 

Theorem 4.1. In an almost unit capacity graph, the D-shortest augmenting path algorithm termi- 
nates in 0{Dm) time. 

The proof of this theorem can be constructed from similar results for the unit capacity case and 
hence omitted. The following theorem provides an approximation bound for the solution produced 
by the /^-shortest augmenting path algorithm. 

Theorem 4.2. Let G{V,E) be an almost unit capacity graph with no parallel arcs. Suppose x is a 
flow generated by the D-shortest augmenting path algorithm and x* be a maximum s — t flow in G. 
Then (i) v{x*) — v{x) < 75^- (H) v{x*) — v(x) < (j^r^)'^- (Hi) If G is a simple almost unit capacity 
graph, then v{x*) — v{x) < 75^2- 
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Again the proof of this theorem can be constructed by modifying arguments in the proof of 
corresponding results for the unit capacity maximum flow algorithms. Detailed proof is available 
in the thesis [63j. 

Let be a flow produced by the D-shortest path algorithm in G. Extend into a maximum 
flow in G using the labeling algorithm [3J. Note that the labeling algorithm performs at most 
v{x*) — v{x^) flow augmentations where x* is a maximum flow in G. Thus this labeling phase can 
be implemented in 0{{v{x*) — v{x^))m) time. Combining this with theorems 14.11 and 14.2( 1) . we get 
a complexity bound of + Dm). Choosing D = \\/rn] , we get a bound of 0(m'^/^). Likewise, 

Combining with theorems 14.11 and I4.2f ii). we get a complexity bound of 0(^^^ + Dm). Choosing 
D = rn2/3], we get a complexity bound of 0{'n?/^m). For simple graphs, we get a complexity 
bound on + Dm). Choosing D = ^/n we get a complexity bound of 0{my/n). Summarizing 

the forgoing discussion, 

Theorem 4.3. The maximum flow problem in almost unit capacity network can be solved in 
0(min{m'^/^, n^/^m}) time. For an almost unit capacity simple graph, the problem can be solved in 
0{my/n) time. 

Note that Theorem 14.31 generalizes corresponding results on unit capacity networks to almost 
unit capacity graphs. The discussions of this section are crucial to our improved algorithm for 
BNFP on unit capacity networks. 

4.2. Algorithm for BNFP in unit capacity networks. Let G{V, E) be a unit capacity graph 
on which a BNFP is deflned. Then, as noted earlier, its auxiliary graph G*{y*,E*) is an almost 
unit capacity graph. Thus combining Theorem 14.31 with the binary search threshold algorithm for 
BNFP discussed in Section 3, it can be seen that BNFP on unit capacity graphs can be solved 
in 0(min{m^/^, n^/^m} logn) time and in 0(m log n-y/n) time on unit capacity simple graphs. We 
now show that we can improve these bounds by a factor of 0{^J\og n). 

Our algorithm first computes a lower bound on the optimal objective function value of BNFP 
using an approximate version of the binary search threshold algorithm. This also generates a flow 
in G* which is a partial solution to BNFP. This partial solution is then extended to a solution to 
BNFP using our BNFP augmenting path algorithm. 

Let Co-(i) < Co-(2) < • • • < Co-((p) be an ascending arrangement of all distinct arc weights of G. 
Note that the optimal objective function value of BNFP is one of these Cj values. Consider the 
graph G'*(Cct(A;)) for some k with edge set, say E'^ . Let x^ be the flow produced by the D-shortest 

augmenting path algorithm on G*(Co-(fc)). If 5 — v{x^) > then by Theorem 14.21 (i), we can 

conclude that the maximum fiow in G*{c^(j.-^) is strictly less than 5 and hence c^(^).-^ is a lower 

bound. If 5 — v{x^) < -j^T^, then the maximum fiow value in G*{c„(^j^^) may or may not be equal 
to 5. In this case, we make a heuristic decision to set Co-(fc) as an approximate upper bound on 
the optimal objective function value. Using this search strategy we present our approximate binary 
search threshold algorithm (Algorithm ABST) below. Without loss of generality assume G* contains 
a maximum fiow of value 6. 

Algorithm ABST 
begin 

construct G* from G 
let c^(i) < c^(2) < 



• • • < Co-(^) be an ascending arrangement of all 
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distinct axe weights of G; 
let Z = 1, u = ip; 
while {u — I > 1) do 

set k:= [^J; 

construct G*{c„(f^^) 

let a;|,_,_2 be the flow produced by the {D + 2)-shortest augmenting 

path algorithm on G*(c^(A;)) = {V*,E*) 
if (5 — v{x'^_^_2) < then u = k; 
else I = k + 1; 
end while; 

(comment: at this stage u = 1.) 

output the flow X produced by the D + 2 shortest augmenting path algorithm on G* {c^^^) ) ■ 



Theorem 4.4. Let x be the flow produced by Algorithm ABST and c = maxjcjj : Xij > 0} = c(o"p). 
Thenc is a lower bound for the optimal objective function value of BNFP. Further, 5 — v{x) < 



Proof. Obviously the starting lower bound Cg-^i) < c*. The index I of the lower bound value 
is updated only when we are guaranteed that G*{c„q^^) contains no flow of value 5 and hence 
C(t(Z+i) ^ c*. Thus by Theorem 3.1, c < c*. Prom the algorithm it can be easily verifled that 

r /-X ^ \E*\ ^ 2\E\ 1-1 



The complexity of Algorithm ABST is O(Dmlogn). If v(x) = S then an optimal solution to 
BNFP can be constructed from x. Otherwise we perform the following extension phase. Starting 
with X and G*(c), the BNPP augmenting path algorithm can be used to compute an optimal 

2 

solution to BNPP by performing at most 0(5) flow augmentations. This can be done in O(^) 
time. Thus BNPP can be solved in 0{Dmlogn + ^) time. Choosing D = ^J^^ we get the 

3 

complexity bound 0(m2 y^Iogn) to solve BNPP on unit capacity graphs. 

A similar algorithm can be obtained if we replace the condition 6 — v{x) < ^-jp hy 5 — v{x) < 
( ' ) in Algorithm ABST. In this case the number of flow augmentation steps in the extension 

/I 2 2 

phase is at most (-^) . This leads to a complexity bound of 0{^^ + Dm log n) time. Choosing 

2 i 2 

D = (j3^) , we get the bound 0(m(n log n) 3 ). The foregoing discussions can be summarized as 

2 3 

Theorem 4.5. The BNFP on a unit capacity graph can be solved in 0(min{m(n log n) s , \/log n}) 
time. 



Note that this bound is better than the bound obtained in the beginning of this section by a 
factor of "v/Iogn. 

Por the case of the bottleneck flows in unit capacity simple graphs, wc replace the condition "if 

l£^*l IV^*I 

5 — v{x) < by "if 6 — v{x) < in algorithm ABST and then as in the previous discussions, 

at the end of Algorithm ABST the extension phase performs at most flow augmentations. Since 
each augmentation takes 0(m) time, it takes O(^) time to complete the extension phase. Thus 

9 



the problem can be solved in + Dm log n) time. Choosing D = y jj^^, we get the complexity 

bound 0{myjn logn). Summarizing these discussions, 

Theorem 4.6. The BNFP on a unit capacity simple graph can he solved in 0{my/n logn) time. 

Note that the bottleneck assignment problem and the bottleneck transportation problem with 
unit capacities are special cases of the BNFP on unit capacity simple networks. Thus we have 
0{m\/n\og n) bounds for these problems as well. This algorithm can be viewed as an extension of 
the algorithm of Gabow and Tarjan [60] for the BAP. 

5. Conclusion 

In this paper we have considered the Bottleneck Network Flow Problem (BNFP), which is a 
generalization of several well studied Bottleneck problems including the Bottleneck Transportation 
Problem (BTP), Bottleneck Assignment Problem (BAP) and Bottleneck Path Problem (BPP). 
Some basic algorithms have been discussed. It is observed that BNFP can be solved as an O(logn) 
sequence of maximum flow problems. Special maximum flow algorithms are identified that can 
easily be modified to solve BNFP in the same worst case time bound as that of solving one maximum 
flow problem by these algorithms. The class of such maximum flow algorithms include generic 
augmenting path algorithms and maximum capacity augmenting path algorithm. We could not 
establish a similar property for the shortest augmenting path algorithm and we have't investigated 
preflow push algorithms in this context. 

We have also considered a special case where the arc capacities are unity. It is well known 

3 2 

that [20] the maximum flow problem on unit capacity graphs can be solved in 0(min{m2 , nsm}) 
time. We showed that the same time bound can be achieved for solving the maximum flow problem 
on an almost unit capacity graph, where the capacities of arcs incident on source and sink nodes are 
allowed to be arbitrary. This together with the binary search threshold algorithms shows that BNFP 

3 2 

on unit capacity networks can be solved in 0(min{m2 , ns m} log n) time. We Also proposed another 

algorithm to solve the problem with an improved complexity of 0(min{m(n log n) s , m 2 ^log n}). 
For the bottleneck flow problems in unit capacity simple graphs, we proposed an 0{m\/n logn) 
algorithm. As a byproduct, we get an 0(n^m logn) algorithm for the bottleneck transportation 
problem with unit capacities. 

An obvious question is if these algorithms can be improved? For the general BNFP, it would be 
interesting to examine what are the maximum flow algorithms that can be modified to solve BNFP 
without increasing the worst case complexity. As noted earlier, BNFP can be solved as a sequence 
of O(logn) maximum flows. Likewise, the maximum flow problem can be solved as an 0(log(nC/)) 
sequence of the BNFPs, where U = maxjujj : {i,j) G E}. It would be interesting to investigate 
further complexity relationships between these problems. 
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